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A COMPUTER PROOF OF A POLYNOMIAL IDENTITY 
IMPLYING A PARTITION THEOREM OF GOLLNITZ 



ALEXANDER BERKOVICH AND AXEL RIESE 



Abstract. In this paper we give a computer proof of a new polynomial iden- 
tity, which extends a recent result of Alladi and the first author. In addition, 
we provide computer proofs for new finite analogs of Jacobi and Euler formu- 
las. All computer proofs are done with the aid of the new computer algebra 
package qMultiSum developed by the second author. qMultiSum implements 
an algorithmic refinement of Wilf and Zeilberger's multi-g-extension of Sister 
Celine's technique utilizing additional ideas of Verbaeten and Wegschaider. 



1. Gollnitz's Partition Theorem and Related 
g-hypergeometric identities 

In 1967, Gollnitz Q proved the following deep partition theorem: 

Theorem 1. Let A(N) denote the number of partitions of N in the form N — 
fii + n-i + 713 + • • • , such that no part is equal to 1 or 3, and ni — Hj+i > 6 with 
strict inequality if ni = 0, 1, 3 (mod 6). 

Let B{N) denote the number of partitions of N into distinct parts = 2, 4, 5 (mod 6). 
Then 

A(N) = B(N). 

In Alladi, Andrews, and Gordon reformulated and refined Theorem |l| using 
the language of colored partitions. To state their theorem, we will need a few 
definitions. 

Let P{N,i,j 1 k) denote the number of partitions of N into parts occurring in 
three (primary) colors ordered as 

A<B<C, (1.1) 

such that parts in the same color are distinct and the number of parts in colors 
A, B, C is equal to i, j, k, respectively. 

Next, consider partitions ir, such that parts equal to 1 may occur in three primary 



colors (1.1), but parts > 2 may occur in six colors ordered as 

AB < AC < A < BC < B < C. (1.2) 

In addition, the gap between adjacent parts is required to be > 1, where the gap 
may equal 1 only if both parts are either of the same primary color or the larger 



part is in a color of higher order according to (1.2). These partitions n were called 
Type-1 partitions in ||. 

We can now state the following result B. 
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Theorem 2. Let G(N, a, b, c, ab, ac, be) denote the number of Type-1 partitions of 
N with a parts in color A, . . . , be parts in color BC. (Note that be is not b ■ c!) 
Then 



E 

constraints 



G(N, a, b, c, ab, ac, be) — P(N, i,j, k), 



where the sum on the left is over the variables a, b, c, ab, ac, be subject to the i,j, k- 
constraints, which here (and throughout) are 



i = a + ab 
j = b + ab 



c ■ 



ac 



ac, 
be, 
-be. 



To see the connection between Theorem |] and Theorem |L we denote part n in 
color A as A„, . . . , part n in color BC as BC„. Next, we replace the colored 
integers A„, . . . , BC„ by regular integers according to the following rules: 



(1.3) 



r A - 


-> 6n — 


4. 




B„ - 


-> 6n — 


2, 




C n - 


-> 6n — 


1, 




AB„ - 


-> 6n — 


6, 


n > 1, 


AC„ - 


-> 6n — 


•5, 


n > 1, 


,BC„ - 


-> 6n — 


3, 


n > 1. 



Note that this replacement converts the ordering 

AB„ < AC„ < A„ < BC„ < B„ < C„, n > 1 

into the conventional ordering 

6n — 6 < 6n — 5 < 6n — 4 < 6n — 3 < 6n — 2 < 6n — 1, n > 1. 

Recalling that part 1 can occur only as Ai,Bi, Ci, we infer that no conventional 
part equals 1 or 3. Also, one can easily check that under the transformations ( |1.3[) 
the color-gap conditions become identical with the gap conditions in Theorem hi 
So, summing over i, j, k one immediately obtains the Gollnitz partition theorem. 

To prove Theorem 0, the authors of || stated it in an analytic form, which they 
called Key Identity, as follows: 



E 



q 



T t +T ab +T ac +T bc . 



(l-q a 



a+bc\ 



constraints 



(q)a(q)b(q)c(q)ab(q)ac(q)bc 



i.q)^i.q)ji.q)k 



(1.4) 



where 



t = a + b + c + ab + ac - 
_ m(m + 1) 



be, 



and for the g-shifted factorial of a is defined as 

fn?=o(i-^), 

(o)„ = (a;g)„ = 



1117=" (i 



aq 



if n > 0, 
if n = 0, 
if n < 0. 
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Several proofs of the Key Identity (1.4) have appeared in the literature [Q, |[ 
[| [tJ. In Alladi and Andrews provided a straightforward g-hypergeometric 
proof of (1.4). Their proof made essential use of Jackson's g-analog of DougalPs 
summation formula. (See, for instance, Gasper and Rahman |5|, (11.21), p. 238].) 
In [Q, Riese used his computer algebra package qMultiSum to find a very simple 
recursive proof of ( |l.4| ). In Q, Alladi and Berkovich proposed and proved a double 
bounded polynomial generalization of (1.4): 

Theorem 3. If i,j,k, L, M are integers, then 

q T t +T ab +T ac +T bc _ 1 

constraints 



L 



L 



-t + a 
a 

-t + a- 
a-1 





~L-t + b~ 




M-t + c 




L-t 




M-t 




M-t 


- <l 


b 




c 


q 


ab 


Q 


ac 




be 



L-t + b 



q s(M+2)-T s +T t - s +T j - s +T k - s 



s>Q 



M-t + c 

c 

L-s 
s,i-s,j - s 





L-t 




M-t 




M-t 


q 


ab 


'! 


ac 


Q 


be - 1 


M - i - 


-J 










k — s 




i 




( 



(1.5) 



In (1.5) and in what follows, the g-binomial and g-multinomial coefficients are 
defined by 



n + m 




n 









0. 



and 



L 




~L~ 




L — a\ 




a%, a,2, . . . , a n _ 


g 




Q 


°2 


9 



if n > 0, 
otherwise, 

L — a\ — a2 — 



— dn-l 



The proof of Theorem ^ given by Alladi and Berkovich J3| used recurrences along 
with Jackson's q-Dougall formula. It is easy to check that in the limit L, M — > oo, 
( pT5| ) reduces to (1.4). Furthermore, if L — M, then the sum on the right in (1.5) 
can be evaluated with the aid of the q-Pfaff-Saalschutz formula. (See, for instance, 
Gasper and Rahman || (11.12), p. 237].) As a result, (1.5) becomes 



E i Tt+Ta ' 



constraints 



x < q 



be 



L 



-t + a 
a 





'L-t + b' 




L-t + c 




L-t 




L-t 




L-t 


- II 


b 


q 


c 


q 


ab 


q 


ac 


q 


be 



L-t 
a 



a-1 
1 





'L-t + b' 




L-t + c 




L-t 




L-t 




L-t 




- '/ 


b 


q 


c 


q 


ab 


q 


ac 


q 


be - 1 


.J 



L - k 




L - i 




L-f 


i 


q 


3 


q 


k 



(1.6) 



A partition theoretical interpretation of (1.6), given in M, can be stated as 
follows: 
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Theorem 4. Let Gl (N, a,b,c, ab, ac, be) denote G(N, a, b, c, ab, ac, be) with the ad- 
ditional constraint that no part exceeds C^. Let P^N,i,j,k) denote P(N,i,j,k) 
with the additional constraints 

A(A)<A L _ fc , A(B)<B L _ 2 , X{C)<C L ^, 

where A(A) is the largest part in color A, and A(B), A(C) have the analogous 



interpretation. Then, for L > max(i + j, j + k,k 

Gl(N 7 a, b, c, ab, ac, be) — Pl(N, i, j, k). 



E 

constraints 



It was pointed out in Q| that (1.6) and Theorem ^ can be employed to derive 
new finite versions of many classical g-hypergeometric identities, including those of 
Gauss, Jacobi and Lebesgue. 

Since Theorem ^] deals with partitions into parts occurring in three primary 
colors, one may suspect that a polynomial anal og of (1.4) should depend on three 
Unitization parameters L\,L 2 ,M and that (L5) is just the special case L\ — L 2 = 
L of this more general result. And, indeed, further i nve stigations led us to the 
following triple bounded polynomial generalization of (|l.4| ) : 



Theorem 5. Let 

g hjtk (Li,L 2 ,M) := 



E 



<i 



Tt+T^+T^+T^! 



constraints 



J"' 



L\ — t + a 





~L 2 -t + b~ 




~L 2 - 1 




~M-t+ c 




~M - 1 




~M - t 


- '1 


b 




ab 


q 


c 


q 


ac 


q 


be 



'Li-t + a-1 




~L 2 -t + b~ 




~L 2 - 1 




~M-t + c 




~M - t 




~M - t 




a- 1 




b 


Q 


ab 


q 


c 


q 


ac 


q 


be - 1 


J 



(1.7) 



and 



Pi,j,k 



{L U L 2 ,M) := 



s>0 



s(M+2)-T s +T i _ a +T 3 _ 3 +T fc _ 



L\ — s 




L 2 - i 




L 2 -i- j + s 




'M-i- f 




i — s 


q 


. J ~ s _ 


q 


s 


q 


k — s 


q 



Then 



g l ,j,k{Li,L 2 ,M) = p i:Jik (Li,L 2 ,M). 



(1.8) 



(1.9) 



We wish t o co mment that if any one of the parameters i,j,k is set to zero, 
then identity ( |l.9| ) reduces to the double bounded key identity for Schur's partition 
theorem B. 



While it is straightforward to extend the analysis of (D to prove (h9), our goal 
here is different. We would like to use (1.9) as the testing ground for the package 
qMultiSum, which has been recently developed by Riese 0. This package is de- 
scribed in some detail in Section 0. In Section [| we will use qMultiSum to obtain 
a nested recursive proof of Theorem |^. In Section || we will give a computer proof 
of the new finite version of Jacobi's formula j^] and then, propose and prove (au- 
tomatically) a new finite analog of Euler's formula. In Section || we will conclude 
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with a brief discussion of the "human insight" , which went into the computer proof 



of (1.9), and with a short description of the prospects for future work. 

2. The Package qMultiSum 

The object of this section is to give a short account on the Mathematica package 
qMultiSum^ which computes recurrences for g-hypergeometric multi-sums. The 
package has been written by the second author and will be described only briefly 
here. For more details the reader is referred to the forthcoming article fij . 

The implementation is based on the method of k-free recurrences, also known 
as Sister Celine's technique (developed by Wilf and Zeilberger |l(J ) . For reasons of 
efficiency we also incorporated ideas from Wegschaider's package MultiSunf] for 
ordinary hypergeometric summation. 

Let n = (m, . . . , n s ) and k = (ki, . . . ,k r ) be vectors of variables ranging over 
the integers. The central concept of (the g-version of) Sister Celine's technique 
is the computation of recurrences for multiple sums ^ k F(n, k), where F(n, k) 
is g-hypergeometric in all of its arguments. For this we proceed by computing a 
so-called k-free recurrence for the summand first. 

Definition 1. A q-hypergeometric function F(n, k) satisfies a k-/ree recurrence, if 
there exist a finite set S of integer tuples of length s + r and polynomials o"ij(n) 
not all zero, such that 

]T a iJ (n)F(n-i,k-j) = (2.1) 

(U)es 



holds at every point (n, k) where all values of F occurring in (2.1) are well-defined. 
The set S is called a structure set. 

The computation of a k-free recurrence is done by making an Ansatz of the form 



(2.1) for some structure set S and undetermined erj.j. Dividing equation (2.1) by 



F(n, k) we get the rational equation 



E 

(ij)es 



<7ij(n)fl Fii j(n,k) =0, (2.2) 



which after clearing denominators turns into the polynomial equation 

J2 oij(n)P J ,, 1J (n,k)=0. (2.3) 
(ij)es 

Next we compare the coefficients of all power products g^ 1 ' 1 • • • q k,Jr in (2.3) with 
zero to get a homogeneous system of linear equations for the <7i.j(n). Note that in 
this system the coefficients are rational functions in several variables and not simply 
numbers. Every non-trivial solution of the system gives rise to a k-free recurrence. 

It has been shown by Wilf and Zeilberger [jl0| that every so-called g-proper hy- 
pergeometric function F(n, k) satisfies a k-free recurrence over some "rectangular" 
structure set Sij := {(i,j) | <«/</;, < j m < J m }, since for large enough I 
and J the number of unknowns exceeds the number of equations. 

However, this result is important only from theoretical point of view, because 
in practice the run-time and memory demand grow very fast with the size of the 



■"-available at http://www.risc.uni-linz.ac.at/research/combinat/risc/software/qMultiSum 
2 available at http://www.risc.uni-linz.ac.at/research/combinat/risc/software/MultiSum 
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structure set S. In particular, it turns out that rectangular structure sets are in 
general not usable, since most recurrences live over a different domain, i.e. many 
(i, j) £ Si j are superfluous points with aij = 0. 

Hence we also generalized the concept of P '-maximal structure sets to the q- 
case, leading to more satisfactory results. The underlying existence theory was 
originally introduced by Verbaeten S for single-sums in the q — 1 case. Since it is 
based on arguments from plane geometry, there is no direct generalization to multi- 
sums. Nevertheless, as Wegschaider (9| pointed out, P-maximal structure sets can 
be computed also in this situation: the idea is to start with a small rectangular 
structure set and then to add all those points that do not increase the degree of the 



polynomial on the left-hand side of (2.3). This way the number of equations in the 
corresponding linear system remains the same, whereas we maximize the number 
of unknowns. 

However, there are still some cases where also this method, called Verbaeten 
completion, misses the minimal structure set. In particular, this happens with 
most of the identities from Section ||. Hence we omit the details (see Riese [Q ) and 
only remark that in these specific instances we could overcome the problem by first 
computing the structure set for q = 1 and using the same for the g-case. Right now 
we do not have an explanation why this actually works. 

We want to emphasize that the correctness of a k-free recurrence computed by 
our program can be checked independently. For that one simply divides equa- 



tion (2.1) by F(n, k) and verifies the resulting rational function identity (2.2) for- 
mally. As Wegschaider showed this implies the correctness of the k-free recurrence 
even at those points (n, k) where F(n, k) = 0. 

Once we have computed a k-free recurrence for the summand, the recurrence for 
the whole sum can be obtained by summing over it. Doing so the left-hand side 
of the recurrence could collapse to 0. Since we do not know of any example where 
this actually occurs, we do not go further into the details here (see Wegschaider || 
or Riese @). 

Moreover, one should keep in mind that for some specific n certain coefficients in 
the recurrence for the sum might vanish. In this case one possibly has to consider 
extra boundaries. 

3. A Computer Proof of the Triple Bounded Identity 



In this section we prove algorithmically identity (l.£) which we restate for con- 
venience as 

gi, j ,k(Li,L 2 ,M)=pi tj , k (L 1 ,L 2 ,M). (3.1) 

For this we proceed in two steps. First we compute a recurrence for ^^(Li, L2, M) 
and show that pi^ y k{L\, L2, M) satisfies the same recurrence. Then we prove that 
( pTl| ) holds at a certain boundary. The latter is achieved by repeating both steps 
for the boundary identity. 

3.1. The Recurrence. First of all we load the package: 

ln[l]:= «qMultiSum.m 
Out[l]= Axel Riese 's qMultiSum implementation version 2.1 loaded 
Then we enter the constraints and the summands. 
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ln[2]:= a = i-ab-ac; b = j-ab-bc; c = k-ac-bc; 
t = a+b+c+ab+ac+bc ; 
T[m_] := m(m+l)/2; 

ln[3]:= gsum = q~(T[t]+T[ab]+T[ac]+T[bc-l]) qBinomial [Li-t+a,a,q] * 
qBinomial [L2 _ t+b,b,q] qBinomial [M-t+c , c , q] * 

qBinomial [L2~t , ab,q] qBinomial [M-t , ac ,q] qBinomial [M-t , be ,q] * 
(q~bc + (l-q-a)/(l-q"(Li-t+a)) (l-q"bc)/ (1-q* (M-t-bc+1) ) ) ; 

ln[4]:= psum = q~(s(M+2)-T[s]+T[i-s]+T[j-s]+T[k-s]) * 

qBinomial [Li-s , i-s ,q] qBinomial [I<2 - i , j-s ,q] * 
qBinomial [L2~i-j+s , s ,q] qBinomial [M-i-j ,k-s ,q] ; 

The function for computing k-free recurrences is called qFindRecurrence (or qFR 
for short). It takes as arguments the summand, the recurrence variables, the sum- 
mation variables, the dimensions of the initial rectangular structure set over which 
Verbaeten completion is performed, and some optional parameters. In our case we 
specify a structure set explicitly, which — as mentioned above — we obtained by 
investigating the q = 1 case. 

In[5]:= qFindRecurrence [gsum, {Li ,L2 ,M, i , j} , {ab, ac.be}, 
{0,0,0,0,0}, {0,0,0}, 

StructSet -> {{0,0,0,0,0,0,0,0}, {0,1,0,0,0,0,0,0}, 
{1,1,1,0,1,0,0,0}, {1,2,1,1,1,1,0,0}}] 

Out[S]= 

q k+L 2 F( -_ 1 + Llj _ 2 + L 2 ,-1 + M,-l + i,-l+ j, -1 + ab, ac, be) + 
q k+L 2 F( ._ 1 +Ll} _ 1 + L ^_ 1 + M ^_ 1+ ^ ab ^ aC] bc j + 

q k F(Li, —1 + 1/2, M, i,j, ab, ac, be) ~ q k F(Li, L 2 ,M, ab, ac, be) = 

Then we sum up this recurrence by applying the function qSumRecurrence (or qSR 
for short) to the previous result. Of course we have to provide the information that 
the first 5 variables are recurrence variables. Since we prefer backward shifts, we 
also call the function BackwardShif ts. 

In[6]:= qSumRecurrence ['/, , 5] // BackwardShif ts 
Out[6]= 

q L2 SXJM(-l + Li,-2 + L 2 ,-l + M,-l+i,-l+j) + 

q L2 SUM(-l + Li,-l + L 2 ,-l + M,i,-l+j) + (3.2) 

SUM{Li,-l + L 2 ,M,i,j) -SVM(Li,L 2 ,M,i,j) =0 

Next we want to show that also Pij,k(Li, Li, M) sa tisfi es the same recurrence. 
Surprisingly, with our package we are not able to find ( |3.2[ ) directly. On the other 
han d ver ification is a trivial task: we simply show that already the summand of p ful- 



fills (3.2) and therefore alsop itself. For this we call the function qCheckRecurrence 
(or qCR for short). 

In[7]:= qCheckRecurrence ['/, , psum] 
Out [7]= True 



3.2. The Boundary. First we note that identity (3J) is true if any of the param- 
eters k is negative, since both sides vanish in this case. Thus, to complete the 
proof it suffices to show that 

g idt k {Li, M) = p itj , k (Li,M), (3.3) 
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where 

gi,j,k(Li,M) gi,j,k(Li,i + j - 1,M) 

and 

p it j t k(Li,M) := Pi,j,k(Li,i+j - 1,M). 

To see t his we have to distinguish two cases. If L 2 > i + j — 1 we rewrite 
recurrence (3.2) as 

S\JM(L u L 2 ,M,i,j) 

= q L2 SUM(-1 + Li, -2 + L a , -1 + M, -1 + 1, -1 + j) + • • ■ 

and find that on the new right-hand side cither j or L 2 is shifted backwards by 1. 
Analogously, if Li < i + j — 1 we rewrite (3.2) as 

SUM(L 1( -l + L 2 ,M,i,j) 

= -q La SUM(-l + L 1 ,-2 + L 2 ,-l + M,-l + i,-l+ j) 

and find that on the new right-hand side cither j is shifted backwards or Li is shifted 
forwards by f. Hence, if we show ( |3.3| ) for all integer valued Li, then identity ( |3.l[ ) 
holds for all integer valued L\ and L 2 . 



To prove (3.3) we proceed as before. This means that we compute a recurrence 
for gij.k(Li, M), check that pi t j t k(Li, M) satisfies the same recurrence and prove 
( f3~3| ) at a certain boundary. 

3.2.1. The Recurrence for the Boundary. This time we succeed without spying out 
the structure set from the q = 1 case: 

ln[8]:= qFR[gsum /. L 2 -> i+j-1, {Xi ,M,i}, {ab.ac.bc}-, -[0,0,0}-, {0,0,1}] // 
qSR[#, 3]& // BackwardShifts 

Out[8]= 

q Ll SUM(-l + L 1 ,-l + Af,-l + i) + SUM(-l + L 1 ,M,i)-SUM(Li,M,i) =0 (3.4) 

ln[9]:= qCR[7,, psum /. L 2 -> i+j-1] 
Out[9]= True 



3.2.2. T he B oundary of the Boundary. Finally, to complete the proof of (3J3) and 
thus of (|3.]|), we prove the identity 

- hM) = Pi,j,k(i + j - (3.5) 

Again we see that this is sufficient by viewing recurrence (|0|) as 

SUM(Li, M, i) = q Ll SUM(-1 + L X) -1 + M, -1 + i) + SUM(-1 + L U M, i) 

in case of L\ > i + j — 1 , and as 

SUM(-1 + Li, M, i) = -q Ll SUM(-1 + Li, -1 + M, -1 + i) + SUM(L X , M, i) 

if L\ < i + j — 1. Hence, if we prove ([3~5|), then identity ( |3 ,3| ) holds for all integer 
valued L\. 

From the single-sum we immediately read off the relation 



Pi,3,k{i + j - 1, M) = 5 ii0 S jfi q 



k 
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where A = M 



j. Hence our boundary identity to verify becomes 
9i,j,k(i +j- 1,M) = S hQ 5j fi q" ^ 



(3.6) 



The reason for switching from M to A here is that proving identity (3.6) for A = 
is easy, whereas for M = it is not at all. 

In[10]:= qFRfgsum /. {Li -> i+j-1, L 2 -> i+j-i, M -> A+i+j}, 
{A,i,j,k}, {ab,ac,bc}, {0,0,0,0}, {0,0,0}, 

StructSet -> {{0,0,0,0,0,0,0}, {0,1,1,1,1,0,0}, {1,0,0,0,0,0,0}, 

{1,0,0,1,0,0,0}, {1,0,1,1,0,0,0}, {1,0,1,1,0,0,1}, 

{1,1,0,1,0,0,0}, {1,1,0,1,0,1,0}, {1,1,1,1,0,0,0}, 
{1,1,1,2,0,1,1}}] // 
qSR[#, 4]& // BackwardShifts 



-2 + ft) 
-1 + fc)- 



(3.7) 



Out[10]= 

9 -l+2A+2 l+ 2i SUM (_ 1 + A, -1 + i, -1 + j 

g- 3+A+3l+3j SUM(-1 + A, -1 + i, -1 + j, 
q A+i+i (-1 + q~ 1+ ' +j ) SUM(-1 + A, -1 + i,j, -1 + ft) - 
q A+ i+j ( _ x + q -i+i+jj SUM (_ 1 + A> i7 _! + j: _i + fc) + 
g A+i+j SUM (_x + A , i,j, -1 + k) + SUM(-1 + A, ft) - 
g -i + A+2 I+ 2 3 SUM (A,-1 +i,-l + +fc) - SUM(A,i,j,fc) = 

Obviously Sjfi q Tk [^] is a solution of this recurrence, we only have to consider 
the two non-degenerate cases when i = j = or i = j = 1. Note that once we have 
proved the validity of (3.6) for A = 0, which follows immediately from [|[ (3.7)], 
our recurrence implies the validity both for A > and A < 0. Indeed, for the case 



A > we rewrite (3.7) as 

SUM(A, k) = g - 1 + 2A + 2 «+ 2 i SUM(-1 + A, -1 + i, -1 + j, -2 + k) + ■ ■ ■ 

and observe that on the new right-hand side at least one of the parameters among 
i, j, k is shifted backwards or A is shifted backwards by 1. Analogously, for the case 
A < we rewrite (3.7) as 

q -l+2A + 2i+2j SUM (_ 1 + 



SUM(-1 + A,i,j,k) 



-l+i,-l+j,-2 + k) 



and observe that on the new right-hand side at least one of the parameters among 
k is shifted backwards or A is shifted forwards by 1. 



Finally, let us briefly summarize the run-times. The computation of (3.7) took 
approximately 35 seconds on an SGI Octane using Mathematica 4.0.1. All other 
computations could be carried out within less than 5 seconds. 



4. Computer Proofs of New Finite Analogs 
of Jacobi and Euler Formulas 

In ||, (5.6)] Alladi and Berkovich pointed out that (1.6) can be used to derive 



-/ 



,22 + 1 



L - k~ 




~L-i 




L - 3 


i 


1 


. j _ 


1 


k 



(4.1) 
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from which Jacobi's formula 

^2(-l) l (2l + l) q T '=( q ;q) 



3 

oo ' 



Z>0 



follows in the limit L — > oo, a — » — 1. It is instructive to prove (4.1) in an automated 
fashion. For the right-hand side of ( |4.l|) we obtain within 20 seconds the following 
recurrence of order 4: 

ln[ll]:= Clear [a] 

ln[12]:= qFR[a~(i-j) (-l)~k q~ (T[i] +T [j] +T [k] ) qBinomial [L-k, i ,q] * 
qBinomial [L-i , j ,q] qBinomial [L-j ,k,q] , 
L, {i,j,k}, 2, {0,0,0}] // qSR 



Out[12] = 



-(ag 9+ ^ SUM(L)) + q +2L (1 - a + a + aq 2+L ) SUM(1 + L) - 
(1 - a + a 2 ) q 4+L (-1 + q 3+L ) SUM(2 + L) + 

(-a - q 4+L + a q 4+L - a 2 q 4+L ) SUM(3 + L) + a SUM(4 + L) = 



Now we plug in the left-hand side of ( 4.1 



ln[13]:= Simplify: 1 /, /. SUM[L + m_] :> SUM[L] + 

Sum[a~(-1) (l+a~(21+l))/(l+a) q~T[l], {l.L+l ,L+m>] ] 

Out[13]= True 



i—3 


L - k~ 




'L-i 




~L-j~ 




i 


q 2 


. j _ 


q 2 


k 



Once we have verified (4.1) for L G {0, 1, 2, 3}, we are done. 
Recently, we came up with the identity 

L 
1=0 

which for L — > oo turns into Euler's formula 

1 = (-9; 9)00 (q;q 2 )oc- 

With our package the proof can be done again in 20 seconds: 

ln[14]:= qFR[(-l)~j q~ (2T [i] +2T [j] +2T [k] -i-j ) qBinomial [L-k, i ,q~2] * 
qBinomial [L-i , j ,q~2] qBinomial [L-j ,k,q~2] , 
L, {i,j,k}, 2, {0,0,0}] // qSR 



Out[14] = 



(q 14+6L SUM(L)) + q 12+4L (1 + q 4+2L ) SUM(1 + L) 



6+2L , 



(-1 + q 3+L ) (1 + q i+1J ) SUM(2 + L) + 



3+L-, 



(-1-9 



&+2L*. 



SUM(3 + L) + SUM(4 + L) = 



Plugging in the left-hand side of (4.2) gives: 

ln[15]:= Simplify ['/. /. SUM[L + m_] : > q~ (2 (L m + T [m] ) ) * 

(SUM[L] + Sum[q~(2(T[L]-T[l] )) , {1 ,L+1 ,L+m}] )] 

Out[lS]= True 



Again the remaining task is to verify (4.2) for L € {0, 1, 2, 3}. 



(4.2) 
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5. Concluding Remarks 



While it is apparent that the package qMultiSum is a valuable tool for proving 
g-hypergeometric identities, we would like to point out that our proof of ( |l.9| ) is not 
fully automated. First of all, in general, the computer generates not just one, but 
many recurrences. Some of these recurrences are "dead ends" in practical terms. 

For instance, if one starts with the simple recursion relation 



g itjfk (L\,L 2 ,M) = gijM{Li - l,L 2 ,M) + q Ll g. l ^ 1 . jtk (L 1 

and then sets L\ = i — 1, one gets 
gi ^ k {i~l,L 2 ,M) = q^ M - 



l,L a -l,M ■ 





L 2 - i 




L 2 - 3 




'M-i- j 




. 3 - * _ 


q 


i 


Q 


k — i 



1) 

(5.1) 



(5.2) 
with 



1.9), tak- 



Unfortunately, we did not succeed in proving the boundary identity (5.2 
qMultiSum and, as a result, we were not able to complete the proof of ( 
ing (5.1) as the starting point. 

Second, crucial "human insight" was used in selecting the not so obvious bound 
aries 

Li=i + j-l, L 2 =i + j-l, M = i+j 



in Section ^2. These boundaries were determined by two requirements: 

1. The validity of the identity on the boundary, along with the recursion rela- 
tions, should imply the validity of the identity everywhere. 

2. The identity should take a particularly simple form on the chosen boundary. 

We have already mentioned at the end of Section || that the recursion relations 
generated by the computer should be examined carefully, because on some specific 
hyperplanes these recurrences may become a triviality = 0. At present, this 
examination has to be carried out manually. 



We would like to finish this article by pointing out that currently only (1.6) 
finds a partition theoretical interpretation. Clearly, it is highly desirable to find a 
partition theorem which corresponds to identity (1.9). It is encouraging that for 



the case when one of the parameters i,j, k in (1.9) is set to zero, such a theorem 
was recently found in . 

Acknowledgment. We would like to thank Krishnaswami Alladi, George E. An- 
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